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We have used domain wall fermions to calculate K —> tt and K ^ matrix elements which can be used to 
study the A7 =1/2 rule for K decays in the Standard Model. Nonlinearities in the A/ = 3/2 matrix elements due 
to chiral logarithms are explored and the subtractions needed for the A7 = 1/2 matrix elements are discussed. 
Using renormalization factors calculated using non-perturbative renormalization then yields values for real 
and A2. We present the details of our quenched 16^ x 32 x 16, /9 = 6.0, M5 = 1.8 simulation, where a previous 
calculation showed that the finite Ls chiral symmetry breaking effects are small (rriros ~ 4MeV). 



1. INTRODUCTION 

At energies below the electroweak scale the 
weak interactions are described by local four- 
fermi operators multiplied by effective coupling 
constants, the Wilson coefficients. The formal 
framework to achieve this is the operator product 
expansion (OPE) which allows one to separate 
the calculation of a physical amplitude into two 
distinct parts: the short distance (perturbative) 
calculation of the Wilson coefficients and the long 
distance (generally non-perturbative) calculation 
of the hadronic matrix elements of the operators 
Qi. We calculate on the lattice K ^ tt and 
K 0. This allows us to calculate the low energy 
constants in chiral perturbation [Q which, after 
incorporating the non-perturbative renormaliza- 
tion factors are then translated into K tttt ma- 
trix elements. The CP-PACS collaboration has 
also presented a very similar calculation at this 
meeting Q. 

2. DETAILS OF THE SIMULATIONS 

We have used the Wilson gauge action, 
quenched, at = 6.0 on a 16'^ x 32 lattice 
which corresponds to an inverse lattice spacing 
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a^^ = 1.922 GeV. The domain wall fcrmion 
height M5 = 1.8 and fifth dimension Lg ~ 
16 give a residual symmetry breaking mrcs = 
0.00124 w 4MeV [|; 400 configurations sep- 
arated by 10000 heat-bath sweeps were used in 
this analysis. K ir matrix elements were cal- 
culated in the S'J7(3)flavor limit for 5 light quark 
masses ruf = 0.01, 0.02, 0.03, 0.04, 0.05. Since the 
K —> matrix elements vanish in the SU{3) limit 
these matrix elements were calculated with non- 
degenerate quark propagators for 10 mass com- 
binations subject to the constraint — rrid = 
0.01, 0.02, 0.03, 0.04. We have also calculated the 
so called eye diagrams with an active charm quark 
for nif = 0.1,0.2,0.3,0.4 ( the physical charm 
quark is around 0.5). However, the analysis for 
charm-in is still in progress; in this presentation 
we concentrate on the case with 3-active flavors 
wherein charm is integrated out assuming it is 
very heavy. The calculation took about 4 months 
on 800 Gflops (peak). Quark propagators were 
calculated using the conjugate gradient method 
with a stopping residual of 10~^ with periodic and 
anti-periodic boundary conditions which amounts 
to doubling the lattice size in time direction. 
The two wall source propagators at Ik = ^ and 
— 27 were fixed to Coulomb gauge. For eye di- 
agrams we employed random wall sources spread 
over time slices t = 14—17 with 2 hits per configu- 
ration. Dividing the three-point correlation func- 
tions by the wall-wall pseudoscalar-pseudoscalar 
correlation function yields the desired matrix ele- 



2 



ments up to a factor of 2m,r which is determined 
from a covariant fit to the wall-point two-point 
function in the range t = 12 — 20 for each mass. 

3. CALCULATION OF LOW ENERGY 
CONSTANTS 

Since our results unambiguously show that 
ReAo and ReA2 come essentially from the 
current-current operators (recall these have the 
largest Wilson coefficients) we will concentrate on 
these operators from now on. Quenched chiral 
perturbation theory predicts 
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We find a quenched chiral logarithm coefficient 
S = 0.05(2) which has a negligible contribution 
in our matrix element calculation. Unlike the 
quenched chiral logarithms, the conventional log- 
arithms coming from quenched chiral perturba- 
tion theory induce large corrections to the A/ = 
3/2 K —>■ TT matrix element as can be seen in 
Figure ||. We fit these amplitudes to jj] 
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'M - f^Trimf + rrires), f = 137 MeV, 
A = 1 GeV. The conventional chiral logarithm 
rnfj \n{m\,^ / A"^) is almost linear over the mass 
range we have used so the fitting routine can- 
not distinguish this term from the linear term 
if we leave the coefficient of the logarithm as a 
free parameter. Since the large coefficient -6 of 
the logarithm makes the contribution of this term 
comparable to the contribution of the linear term 

(27 1) 

omitting this term would change bl by almost 
a factor of two. The quenched chiral log contri- 
bution is very small. A/ = 1/2 K tt matrix 
elements mix with sd with a power divergent co- 
efficient ^ {nis + rad)/a? ■ We define a subtracted 
matrix element {'K^\Q'f''^^ \K^) sub by 

{Ti+\Qf'^'^\K+) + 77,(m, + md){'K+\M\K+) 

where 77, is obtained from a linear fit to . 
For an explanation of this subtraction of the 
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Figure 1. The matrix element (7r+|e|f J'^^ 

power divergence we refer the reader to The 
quenched chiral perturbation theory corrections 
to (tt+IQ^^I are not known. As can be 

seen in Figure ^ our data is consistent with a lin- 
ear fit co^i + Ci^iiTif with the slope Ci^i determining 
the low energy constants and the intercept 
co,i arising from residual chiral symmetry break- 
ing. 
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Fiffure 2. The matrix element (n+\Q/ \K+)sub 



4. Re ^0 AND Re A2 

We use chiral perturbation theory to compute 
the lattice K — > tttt matrix elements. Using 
non-perturbative Z factors we obtain the contin- 
uum matrix elements which are then multiplied 
by Wilson coefficients to yield the physical am- 
plitudes. We present an extrapolation to the 
Kaon mass scale to lowest order in chiral pertur- 
bation theory and a second extrapolation which 
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includes one loop logarithmic effects. We multi- 
ply the pseudoscalar masses by ^ so that for ^ — > 
the chiral perturbation theory extrapolation is in- 
creasingly accurate but we need the extrapolation 
at ^ = 1, the physical point. In figure^ we present 
HeAo and ReA2 as a function of the parameter ^. 
The chiral logarithm correction for Re^o is large 
(about 40%). In addition one expects a large 
correction (not included here) coming from the 
tree level 0(p^) terms necessary to cancel the de- 
pendence on the chiral perturbation theory scale 
A^PT- If the Z factors and the Wilson coefR- 
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Figure 3. KeAo and KeA2 as a function of ^ 

cients were calculated to all orders in perturba- 
tion theory the physical amplitudes that we cal- 
culate would not depend on the scale fi where the 
transition between the lattice and the continuum 
operators is made. To a good approximation this 
is what we find, even though at ^ = 1.51 GeV 
one expects non-perturbative effects in the Z fac- 
tors and at fj, = 3.02 GeV the discretization errors 
may be large. In Figure ^ we present the ratio 
Reylo/Re^2, the so called A/ = 1/2 rule which 
shows a large enhancement in the / = channel 
in accord with experiment (note, the chiral log- 
arithm corrections largely cancel in the ratio so 
a large enhancement is seen for both extrapola- 
tion choices). The residual scale dependence in 
the physical amplitudes is slight (see Figure 0) . 
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Figure 4. RcAq and KeA2 obtained using 1-loop 
logarithmic corrections vs. the matching scale n; 
the dotted lines show the experimental results. 

5. CONCLUSIONS 

In conclusion, KbAq, Re^2 and especially the 
ratio ReAo/ReA2 were found reasonably close to 
the experimental values. We see this as an im- 
portant success of the lattice method. However 
there were a number of major approximations in 
our calculation, the hardest to quantify is the use 
of quenched QCD. Also the chiral logarithms in 
quenched if — > tt. A/ = 1/2 are not known and 
we have included only the logarithmic portion 
of the next-to-leading-order, 1-loop corrections in 
K TTvr extrapolations. 
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